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Abstract 

Background: Keteleeria davidiana var. formosana (Pinaceae), Taiwan cow-tail fir, is an endangered species listed on 
the IUCN Red List of Threatened Species and only two populations remain, both on the Taiwan Island. Sixteen 
polymorphic microsatellite loci were developed in an endangered and endemic gymnosperm species, Keteleeria 
davidiana var. formosana, and were tested in an additional 6 taxa, K. davidiana var. calcarea, K. davidiana var. 
chienpeii, K. evelyniana, K. fortune!, K. fortune! var. cyclolepis, and K. pubescens, to evaluate the genetic variation 
available for conservation management and to reconstruct the phylogeographic patterns of this ancient lineage. 

Findings: Polymorphic primer sets were developed from K. davidiana var. formosana using the modified AFLP and 
magnetic bead enrichment method. The number of alleles ranged from 3 to 1 6, with the observed heterozygosity 
ranging from 0.28 to 1.00. All of the loci were found to be interspecifically amplifiable. 

Conclusions: These polymorphic and transferable loci will be potentially useful for future studies that will focus on 
identifying distinct evolutionary units within species and establishing the phylogeographic patterns and the 
process of speciation among closely related species. 
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Findings 

Background 

For conservation management, the evaluation of genetic 
diversity and the identification of distinct evolutionary 
units in endangered taxa should allow for the develop- 
ment of more efficient conservation strategies [1-3]. 
Over the past century, human activities associated with 
the lumber industry and agriculture have destroyed and 
fragmented forest habitats at altitudes of 500-2500 m 
on Taiwan Island. Several gymnosperm species have 
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decreased in population size and are now under threat, 
including Cycas taitungensis [4-6], Podocarpus nakaii 
[1], Amentotaxus formosana [2], and Keteleeria davidi- 
ana var. formosana [7]. Extreme population declines 
can affect the genetic diversity of a species [8] . The loss 
of genetic diversity carries serious evolutionary con- 
cerns regarding both the ability to adapt to changes on 
oceanic islands and longer-term interactions with other 
organisms [9]. Thus, molecular markers are used to 
evaluate endangered species in greater depth, with the 
purpose of applying this genetic diversity information to 
the conservation and restoration of biodiversity [10]. 

Keteleeria davidiana var. formosana (Pinaceae), Taiwan 
cow-tail fir, is an endangered species listed on the IUCN 
Red List of Threatened Species [11] and derived from K. 
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davidiana, including var. davidiana and var. calcarea, 
occurred in China by geographical disjunction and leaf 
scars obscurely or obviously protruding on branchlets 
[12]. The taxa has only two wild populations, the Dawu 
Taiwan Keteleeria Forest Reserve (Dawu) and the Pin- 
glin Taiwan Keteleeria Nature Reserve (Pinglin), which 
are located in Northern and Southern Taiwan, respect- 
ively [13]. A previous inventory study has found that 
less than a hundred individuals remain in the wild, few 
seedlings are found in the field, and adult trees are vulner- 
able to lightning damage or death [14]. Additionally, con- 
tinuous logging over the past century in the low-elevation 
broadleaf forests has restricted K. davidiana var. formosana 
to narrow habitats and has decreased its population size 
[15]. Therefore, the conservation strategies for Keteleeria 
davidiana var. formosana should include in situ protection 
and management in conjunction with ex situ approaches 
such as clonal orchards. Hence, the population genetics of 
this species should be evaluated before in situ and ex situ 
conservation management begins. In this study, we devel- 
oped polymorphic microsatellite loci from Taiwan cow-tail 
fir to estimate the population structure within the species 
and to identify distinct evolutionary units within the popu- 
lation. In addition, the transferability of the microsatellite 
primers developed in this study was tested in other taxa, 
specifically, three species and three varieties of the Kete- 
leeria genus. 

Methods 

Sampling and DNA extractions 

Twenty-five individuals were sampled from each of the 
two remaining populations of K. davidiana var. formo- 
sana from Dawu and Pinglin, Taiwan (Table 1). To test 
the transferability of the markers, two to three individ- 
uals of three species, K. pubescens, K. fortune, and K. 
evelyniana, and three varieties, K. fortunei var. cyclolepis, 
var. calcarea, and K. davidiana var. chienpeii, were sam- 
pled from the field. Permission for tissue collection was 
obtained from the Council of Agriculture, Republic of 
China, authorities. All taxa of the Keteleeria genus are 



diploid with 24 chromosomes [13]. The sample size, lo- 
cation, and storage herbarium for the voucher specimens 
are listed in Table 1. Total genomic DNA was extracted 
from silica-dried leaf powder using the Plant Genomic 
DNA Extraction Kit (RBC Bioscience, Taipei, Taiwan). 

Isolation of microsatellite DNA loci and identification 

To isolate the microsatellites efficiently, we used the modi- 
fied AFLP [16] and magnetic bead enrichment method 
[1,17,18] to select microsatellite loci. Genomic DNA from 
one individual of K. davidiana var. formosana was treated 
with the restriction enzyme Msel (Promega, Madison, 
Wisconsin, USA) and separated using 1% Nusieve" 3:1 
agarose gel (FMC Bio Products, Rockland, ME, USA) 
electrophoresis. The cut DNA fragments ranged from 
400 to 1000 bps, and the purified fragments were 
extracted from the agarose gel using the HiYield™ Gel 
PCR DNA Fragments Extraction Kit (RBC Bioscience). 
Purified DNA fragments were ligated to a double-stranded 
Msel-adaptor (complementary oligo A: 5 -TACTCAG 
GACTCAT-3'; 5 phosphorylated oligo B: 5'-GACGAT 
GAGTCCTGAG-3 ') using the Quick Ligation™ Kit 
(New England Biolabs, MA, USA) at 25°C for 5 minutes. 
The ligation product was used as the template DNA for the 
enrichment of the partial genomic library, and the adaptor- 
specific primers, named Mse I-N (5-GATGAGTCCT 
GAGTAAN-3'), were used to perform 20 cycles of prehy- 
bridization PCR amplification. The PCR mixture contained 
20 ng template DNA, 10 pmol adapter-specific primer, 
2 uL 10 x reaction buffer, 2 mM dNTP mix, 2 mM MgCl 2 , 
0.5 U Taq DNA polymerase (Promega), and sterile water 
was added to reach a total volume of 20 uL. The PCR 
protocol was set at 94°C for 5 min, followed by 18 cycles of 
94°C for 30 s, 53°C for 1 min, and 72°C for 1 min, using a 
Labnet MultiGene 96-well Gradient Thermal Cycler 
(Labnet, Edison, New Jersey, USA). The amplicons 
were denatured and hybridized to four separated 
5 -biotinylated oligonucleotide probes, (AG)is, (AC)is, 
(TCC)io, and (TTG)i 0 , in 250 uL of hybridization 
solution at 68°C for 1 h, and 1 mg of Streptavidin 



Table 1 Sample location for each species, variety, and population of the Keteleeria species 



Species 


Location 


Sample size 


Latitude 


Longitude 


Herbarium 


K. davidiana var. formosana 


Pinglin, Taiwan 


25 


N 24°53'28" 


E 


121°42'12" 


TAIF 




Dawu, Taiwan 


25 


N 22°2V23" 


E 


120°47'27" 


TAIF 


K. fortunei var. cyclolepis 


Fujian, China 


3 


N 27°09'18" 


E 


118°13'03" 


TAIF 


K. pubescens 


Guizhou, China 


3 


N 26°30'54" 


E 


108°03'52" 


TAIE 


K. davidiana var. calcarea 


Guangxi, China 


3 


N 25°12'47" 


E 


1 1 0°1 1 '50" 


TAIE 


K. davidiana var. chienpeii 


Guizhou, China 


3 


N 26°20'26" 


E 


106°41"52" 


TAIE 


K. fortunei 


Fujian, China 


3 


N 26°09'57" 


E 


1 19°02'19" 


TAIF 


K. evelyniana 


Yunnan, China 


3 


N 23°42'55" 


E 


105 o 07'24" 


TAIE 



Sample size, location, coordinates, and voucher specimens are indicated. 

Note: TAIF = the herbarium of the Taiwan Forestry Research Institute; TAIE = the herbarium of the Taiwan Endemic Species Research Institute. 
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MagneSphere Paramagnetic Particles (Promega) was 
added to the mixture solution to capture the hybridiza- 
tions at 42°C for 2 h. The four enriched microsatellite 
DNA fragments were washed with high-salt and low-salt 
solutions and then used as templates for 25 cycles of PCR 
amplification using the adaptor-specific primers; the ampli- 
fication protocol was same as that used for the prehybridi- 
zation PCR. The PCR products were purified using the 
HiYield™ Gel PCR DNA Fragments Extraction Kit (RBC 
Bioscience) and then cloned using the pGEM 8 -T Easy 
Vector System (Promega). The plasmid DNAs that were 
purified from white colonies were digested by the restric- 
tion enzyme £coRI (Promega) and screened using 1% agar- 
ose electrophoresis to determine the size of the inserted 
DNA fragments. In total, 670 white colonies were selected 
to purify plasmids for sequencing from Keteleeria davidi- 
ana var. formosana. The sequencing cocktail contained 
200 ng plasmid DNA as the template DNA, 0.5 mM T7 or 
SP6 primers, 4 m.L Big Dye Terminator v3.1 Ready Reaction 
Premix Reactions (Applied Biosystems, Carlsbad, Califor- 
nia, USA), and sequencing buffer in a total volume of 
20 uL. Selected plasmid DNAs were sequenced in both 
directions using an ABI PRISM 8 3700 DNA Sequencer 
(Applied Biosystems). These sequences, including their 
microsatellite loci, were recognized using Tandem Re- 
peats Finder version 4.07b [19], and a pair of specific 
primers for each microsatellite locus was designed using 
FastPCR software version 6.4.18 [20]. 

DNA amplification and genotyping 

The optimal annealing temperature was evaluated by gradi- 
ent PCR using a Labnet MultiGene 96-well Gradient Ther- 
mal Cycler (Labnet) for a temperature range from 50°C to 
65°C. For each Keteleeria taxon, two to three individuals 
were evaluated for the optimal annealing temperature. The 
polymorphism of the two remaining populations of K. 
davidiana var. formosana was evaluated using 25 individ- 
uals from each population (Table 1). To test the optimal an- 
nealing temperature, a 20 uL reaction mixture containing 
20 ng template DNA, 0.2 \iM each of the reverse and for- 
ward primers, 2 \iL 10 x reaction buffer, 2 mM dNTP mix, 
2 mM MgCl 2 , 0.5 U Taq DNA polymerase (Promega), and 
sterile water was amplified using a Labnet MultiGene 96- 
well Gradient Thermal Cycler (Labnet). The PCR program 
was set at 94°C for 5 min, followed by 30 cycles of 94°C 
for 40 s, a temperature gradient ranging from 50 to 
65°C for 60 s, and 72°C for 60 s, and a final extension of 
72°C for 10 minutes [21]. Amplicons were checked by 
1% agarose electrophoresis to isolate the target DNA 
bands, which were confirmed by sequencing [18]. To 
examine the genetic polymorphisms, the PCR cocktail 
was amplified using the previously described program 
with the temperature gradient replaced by the optimal an- 
nealing temperature (Ta). The amplicons were checked by 



1% agarose electrophoresis to isolate the target DNA 
bands, which were confirmed by sequencing [18]. The 
amplicons were separated by electrophoresis on a 10% 
polyacrylamide gel (acrylamide: bisacrylamide 29: 1, 80 V 
for 14-16 hours) using a 25 or 50 bp DNA Step Ladder 
(Promega) to determine the allele size. The bands were 
then imaged under UV light using the Flo Gel FGIS-3 
fluorescent gel image system (Top BIO Co., Taipei, Taiwan), 
and the sizes of the PCR products were identified using 
Quantity One software version 4.62 (Bio-Rad Laboratories, 
Hercules, California, USA). 

Data analysis 

The genetic variation indices, including the number of al- 
leles (A/a), the number of effective alleles (Ne), the ob- 
served heterozygosity (Ho) and expected heterozygosity 
(He), Shannon's information index (H), and fixation index 
(F IS ) were estimated using GenAlEx version 6.5 [22]. The 
Hardy- Weinberg equilibrium (H WE ) was tested using 
Arlequin software version 3.5.1.2 [23]. 

Results and discussion 

Development of polymorphic microsatellite markers 

A total of 392 microsatellite loci containing repeat motifs 
were discovered, with maximum and minimum lengths of 
992 bps and 119 bps, respectively, and an average sequence 
length of 477 bps. We designed 102 primer pairs between 
the up- and down-flanking regions of the motifs based on 
the primer design parameters computed using FastPCR 
software version 6.4.18 [20]. To test the optimal an- 
nealing temperatures, which were obtained using gra- 
dient temperature PCRs, template DNA was derived 
from three individuals of each of the Keteleeria taxa, 
which are listed in Table 1. We tested the polymor- 
phisms for K. davidiana var. formosana using 25 indi- 
viduals from each population. Finally, we selected 16 
polymorphic loci from the 102 microsatellites based on 
the detection of unambiguous polymorphic amplicons 
using fixed annealing temperature PCR and a poly- 
acrylamide gel genotyping protocol. The characteris- 
tics of the 16 polymorphic microsatellite loci are listed 
in Table 2. Of the 16 loci, 11 are complete microsatel- 
lite loci, consisting of 5 with a dinucleotide motif, 5 
with a trinucleotide motif, and 1 with a tetranucleotide 
motif. Of the 5 remaining loci, 3 carried a compound 
motif and 2 carried an interrupted motif. The sequences 
of 16 loci reported in this paper are available from GenBank 
(accession numbers: HG518488-HG518503) (Table 2). 

Genotyping and population genetics analysis 

To examine the extent of genetic polymorphism at each 
locus, 25 individuals were collected in the field from 
each of the two remaining populations of K. davidiana 
var. formosana (Table 1). For the 16 new polymorphic 
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Table 2 Summary of general information for the 16 polymorphic microsatellite loci isolated from Keteleeria davidiana 
var. formosana and their transferability to six Keteleeria taxa 



Locus 


Repeat motif 


Primer sequence(5'-3') 


Allele 


Ta(°Q 














Genbank 


Cross 










size 
(bps) 


TW 




pu 


ca 


CP 


fo 


ev 


accession 
no. 


amplification 


TW-5 




F: 
R: 


TTGTCTAGGTATTTGTGCCC 
AAACGCAGTATAAGACAGGC 


116 


54 


54 


54 


54 


54 


54 


54 


HG5 18488 


7/7 


TW-49 


(CTCC) 9 


F: 
R: 


TCACTCCCTATATCCCTAGC 
GGAGAGAGAGATTGTGTTTG 


99 


50 


50 


50 


50 


50 


50 


50 


HG5 18489 


7/7 


TW-69 


(CT) 27 (TATTTGA) 3 F: 


CCTATATCCTCCCCCACCT 


260 


54 


54 


5---1 


54 


54 


54 


54 


HG5 18490 


7/7 






R: 


AGGTATAGAGGGGGATAGTG 






















TW-107 


(CA) 4 (AG), 4 


F: 
R: 


AACCAAGGATGAAACCCTAG 
AAGTGAAAGACATAGTGAGG 


153 


54 


54 


54 


54 


54 


54 


54 


HG5 18491 


7/7 


TW-172 


(GA) 30 


F: 
R: 


AGACAAAGGAGAATGGTGGG 
CCATGCCATGGGTAAGATAG 


165 


50 


50 


50 


50 


50 


50 


50 


HG5 18492 


/// 


TW-195 


(CT) 18 


F: 
R: 


ATCTTCATCACTCTCCCACC 
AAGGGGAAAGTGAGAGAGGG 


170 


54 


54 


54 


54 


54 


54 


54 


HG5 18493 


/// 


KE-6 


(AC) 15 


F: 
R: 


ACGCGAGGTTTGGGGTGGGA 
CTGGGTACTGGTGCCCTGT 


192 


54 


54 


54 


54 


54 


54 


54 


HG5 18494 


/// 


KE-54 


(AC) 15 


F: 
R: 


ACCTTGCAAAGTCACCAGCT 
CTATGGCTAGGTTTCCTGGGT 


271 


52 


52 


52 


52 


52 


52 


52 


HG5 18495 


/// 


TCC-6 


(AAC) 9 


F: 
R: 


TGTATGGTGTGTCAGATGC 
CCTCCTTGTGGAGTTGGCA 


342 


56 


56 


56 


56 


56 


56 


56 


HG5 18496 


111 


TCC-31 


0TG) 5 0TC) 6 


F: 
R: 


CTTCAAAGCTGCTCGACTC 
GTTTGCTTGAGGAACAGGT 


285 


55 


55 


55 


55 


55 


55 


55 


HG5 18497 


111 


TCC-55 


(GAA) 6 (CAA) 12 


F: 
R: 


CTCTTGAAGAACAGGTTGAG 
CTTCAAAGCTGCTCGGGTGT 


288 


54 


54 


54 


54 


54 


54 


54 


HG5 18498 


111 


TCC-139 


(GAA) 13 


F: 
R: 


CCTACAAGGTTGACCCCGT 
TGCTCGCCAATGCATGACAAG 


220 


54 


54 


54 


54 


54 


54 


54 


HG5 18499 


III 


TCC-245 


(GAA) 13 


F: 
R: 


AGGTrrGAGGAGCAGGTTG 
TGAGAGATGTTCCTG1TCTG 


255 


50 


50 


50 


50 


50 


50 


50 


HG5 18500 


in 


TTG-11 


(CAA), 0 


F: 
R: 


GIN I G I AGCTCTGTTGGAC 
CCATATGACCTAGCTTCCCA 


251 


50 


50 


50 


50 


50 


50 


50 


HG518501 


in 


TTG-49 


(AGG) 8 


F: 
R: 


GAGATGACAGAGGCTGCTG 
GCTAGTGGGTGCACTAGGT 


224 


55 


55 


55 


55 


55 


55 


55 


HG5 18502 


in 


TTG-54 


(GA) 15 (GAA) 2 


F: 
R: 


CTCGTGCTAACCCTTGATC 
GGCTACCTTAGGATGTGTGCA 


222 


50 


50 


50 


50 


50 


50 


50 


HG5 18503 


in 



Note: TW = K. davidiana var. formosana; cy = K. fortunei var. cyclolepis; pu = K. pubescens; ca = K. davidiana var. calcarea; cp = K. davidiana var. chien-peii; 
fo = K. fortune; ev = K. evelyniana; F = the forward primer; R = the reverse primer; Ta = optimized annealing temperature. 



microsatellite loci, the number of alleles per locus (Na) 
and the number of effective alleles (Ne) ranged from 5 
to 14 and 2.59 to 11.25 in the Dawu population, respect- 
ively, and from 3 to 16 and 2.14 to 10.59 in the Pinglin 
population, respectively (Table 3). The observed and ex- 
pected heterozygosity {Ho and He, respectively) varied 
from 0.44 to 1.00 (average of 0.68) and 0.61 to 0.91 
(average of 0.82) in the Dawu population, respectively, 
and from 0.28 (average of 0.63) to 0.80 (average of 0.78) 



in the Pinglin population. The means of Shannon's infor- 
mation index (H) and the fixation index (F IS ) were 1.93 
and 0.175 and 1.73 and 0.188 in the Dawu and Pinglin 
populations, respectively (Table 3). Significant deviations 
from Hardy- Weinberg equilibrium {H WE ) were detected 
at 6 and 3 loci in the two populations, respectively, and 
these deviations were attributed to the heterozygote defi- 
ciency of the endangered species (Table 3). A total of 52 
and 29 private alleles were observed in the Dawu and 
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Table 3 Genetic diversity characteristics of the 16 polymorphic microsatellite loci tested on two populations of 
Keteleeria davidiana var. formosana 

K. davidiana var. formosana 



Dawu population (N = 25) Pinglin population (N = 25) 



Locus 


A/a 


We 


Ho 


He 


H 


F,s 


Wo 


We 


Ho 


He 


H 


F,s 


TW-5 


5 


2.59 


0.44 


0.61 


1.18 


0.284 


3 


2.68 


0.56 


0.63 


1.03 


0.106 


TW-49 


/ 


3.93 


0.68 


0.75 


1.55 


0.088 


3 


2.14 


0.48 


0.53 


0.90 


0.100 


TW-69 


8 


5.30 


0.68 


0.81 


1.84 


0.162 


/ 


4.73 


0.64 


0.79 


1.71 


0.189 


TW-107 


/ 


4.36 


0.56 


0.77* 


1.68 


0.273 


4 


3.24 


0.52 


0.69 


1.28 


0.248 


TW-172 


6 


4.63 


0.68 


0.78 


1.63 


0.133 


3 


2.35 


0.52 


0.57 


0.97 


0.093 


TW-195 


5 


4.96 


0.72 


0.80 


1.61 


0.098 


6 


5.00 


0.72 


0.80 


1.66 


0.100 


KE-6 


12 


9.69 


0.60 


0.90* 


2.37 


0.331 


8 


7.27 


0.28 


0.86* 


2.03 


0.675 


KE-54 


10 


7.62 


0.68 


0.87 


2.14 


0.217 


8 


7.23 


0.68 


0.86 


2.02 


0.211 


TCC-6 


10 


5.79 


0.80 


0.83 


1.98 


0.033 


10 


8.56 


0.80 


0.88 


2.22 


0.094 


TCC-31 


8 


5.76 


0.72 


0.83 


1.89 


0.129 


9 


5.73 


0.72 


0.83 


1.96 


0.128 


TCC-55 


10 


8.50 


0.68 


0.88 


2.22 


0.229 


10 


8.22 


0.76 


0.88 


2.18 


0.135 


TCC-139 


6 


4.92 


0.64 


0.80* 


1.67 


0.197 


10 


5.73 


0.56 


0.83* 


1.96 


0.322 


TCC-245 


8 


6.04 


0.72 


0.83 


1.91 


0.137 


4 


3.33 


0.72 


0.70 


1.29 


-0.029 


TTG-11 


13 


7.76 


0.76 


0.87* 


2.28 


0.128 


8 


6.19 


0.72 


0.84 


1.93 


0.141 


TTG-49 


14 


11.26 


1.00 


0.91* 


2.51 


-0.097 


10 


6.87 


0.80 


0.85 


2.06 


0.064 


TTG-54 


13 


9.92 


0.48 


0.90* 


2.42 


0.466 


16 


10.59 


0.52 


0.91* 


2.55 


0.426 


Mean 


8.88 


6.44 


0.68 


0.82 


1.93 


0.175 


7.44 


5.62 


0.63 


0.78 


1.73 


0.188 



*Significant deviation from Hardy-Weinberg equilibrium: P<0.05. 

The number of different alleles (Na), number of effective alleles [Ne), observed heterozygosity (H 0 ), expected heterozygosity (He), Shannon's information index (H), 
and fixation index (F /5 ) are reported. 



Table 4 Result of cross-species transferability in six Keteleeria taxa using the 16 microsatellite primers developed from 
Keteleeria davidiana var. formosana 

Locus K. fortunei var. K. pubescens K. davidiana var. K. davidiana var. K. fortunei (N = 3) K. evelyniana (N = 3) 

cyclolepis (N = 3) (N = 3) calcarea (N = 3) chienpeii (N = 2) 



TW-5 


I 


2 


1 


2 




1 


TW-49 


2 


2 


2 


1 




2 


TW-69 


3 


2 


2 


2 




2 


TW-107 


3 


1 


2 


2 




1 


TW-172 


2 


1 


2 


1 




2 


TW-195 


2 


2 


3 


2 


2 


1 


KE-6 


2 


3 


3 


2 


2 


2 


KE-54 


3 


3 


2 


1 


2 


2 


TCC-6 


3 


3 


3 


2 


3 


3 


TCC-31 


2 


3 


2 


2 


2 


1 


TCC-55 


3 


2 


1 


2 


1 


2 


TCC-139 


2 


3 


2 


1 


2 


2 


TCC-245 


2 


2 


2 


2 


2 


1 


TTG-11 


3 


3 


3 


2 


3 


3 


TTG-49 


3 


2 


3 


2 


3 


3 


TTG-54 


2 


2 


3 


2 


2 


2 



For loci that were successfully amplified, the number of alleles is given. 
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Pinglin populations, respectively, revealing population 
differentiation between the two remaining populations 
of the species. 

The genetic variability, including the means of the ob- 
served and expected heterozygosity (Table 3), of K. 
davidiana var. formosana was high compared with that 
of other Pinaceae species such as Pinus koraiensis (0.38 and 
0.57) [24] and P. massoniana (0.27 and 0.65) [25,26], but 
was similar to that of P. pinaster (0.65 and 0.83) [27]. Un- 
fortunately, no data for other Keteleeria taxa are available 
for the comparison of genetic variability. However, the high 
observed and expected heterozygosity values implied that 
the species was historically larger in population size and 
was more common in low-elevation forests than it is at 
present [7]. 

Test of transferability 

To test the transferability of these microsatellite loci, we 
tested these primers in six other cow-tail fir taxa: K.fortunei 
var. cyclolepis, K. pubescens, K. davidiana var. calcarea, 
K. davidiana var. chienpeii, K. fortune, and K. evelyniana 
(Table 1). Two to three samples of each taxon were used in 
the evaluation of cross-amplification. All loci were transfer- 
able to the six taxa of Keteleeria (Table 4), and the anneal- 
ing temperatures are listed on Table 2. However, few alleles 
were recovered in the six other Keteleeria taxa sampled 
here were caused by a small sampling effect. Nonetheless, 
the systematic and phylogeographic patterns among the 
Keteleeria species remain unclear. The transferability of 
these microsatellite loci among different Keteleeria taxa in- 
dicates the usefulness of these molecular genetic markers 
for interspecific and intraspecific research on such topics as 
phylogeography, speciation and introgression. 

Conclusions 

For conservation purposes, 16 new polymorphic microsat- 
ellite loci were developed from K. davidiana var. formo- 
sana. The genetic variation indices evaluated using these 
16 polymorphic microsatellite loci for the two remaining 
populations of this endangered species are potentially use- 
ful for future studies that will focus on identifying distinct 
evolutionary units within the populations for conservation 
management. The interspecies transferability of these 
microsatellite loci may also be useful for future research 
aiming to reconstruct the phylogeographic patterns and the 
process of speciation among closely related species. 

Availability of supporting data 

The microsatellite sequences are available through the 
National Centre for Biotechnology Information (see 
http://www.ncbi.nlm.nih.gov/). The accession numbers 
on the repository are the following: GenBank accession 
number HG518488 through HG518503. 
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